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ON DIFFERENTIAL EQUATIONS WITH PERIODIC INTEGRALS. 
By Dr. G. W. Hill, Washington, D. C. 

The independent variable being conceived as time, a sy.stem of differential 
equations may be said to admit periodic integrals when the values of the depen- 
dent variables either exactly, or with approximate tendency, after a certain lapse 
of time, repeat their series of values. In the latter case the larger the lapse is 
made the more nearly is the repetition brought about. Strange as it may seem, 
this subject, except in the case of simply periodic integrals, is, at present, not 
completely understood. The text-books on differential equations are almost 
wholly engaged with the cases in which, by certain artifices, the integration can 
be accomplished in finite terms or reduced to quadratures. In the treatment of 
physical problems, however, equations of this sort are rarely met with. Far 
more frequently it is found that methods of approximation must be resorted to. 
Cauchy appears to be the author who has done most for the elucidation of this 
part of the subject. His memoirs are in his later Exercises and in the volumes 
of the Coniptes Rendus for 1856 and 1857. In this article I propose to show how 
simply periodic integrals arise and afterwards to illustrate the general theory by 
treating a problem relating to the motion of a system of points. 

I. 
Having the independent variable /, and the two dependent variables x and .r, 
let us suppose the latter satisfy the equations 

dx dx, ^, , 

A cross multiplication between the members of these equations gives 

dx, ., .dx 

The integral of this is, C being the arbitrary constant, 

x^^==2ji{x)dx-\- C. 

The values of x and x^ being known for a given value of/, we readily find the 
value of C proper to the special case we treat. By substituting the value of x 
derived from this equation in the first of the differential equations we get 

dx 
di 



■.^\2Ji{x)dx^c\. 
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The expression under the radical sign is a function of x; calling it X, let us con- 
sider the equation X = o. Since real values of x are supposed to correspond 
to all values of t, X can never be negative ; and, from the way the constant C was 
determined, it is plain that, for the given value of t, X is positive. Then in X= o, 
let X be supposed to increase until a value x=^b\?, reached for which X= o, 
that is to say a real root of this equation. Similarly let x diminish from the same 
point until a value ;ir = ^ is reached for which again X^o, that is a second real 
root Then, X being positive for all values of x which lie between c and b, if the 
latter are non-multiple roots, X\% negative for values of x which lie just outside 
these limits. Thus x must necessarily remdin within the limits c and b. Also, 
in its motion, it always attains them; for suppose x is augmenting, then the rad- 
ical, which forms the value of dxjdt, must be taken positively, and, from the law 
of continuity, must continue to be so taken until it becomes zero, that is until x ar- 
rives at the value b. But dxjdt cannot be positive beyond this point, for x can- 
not surpass b. Hence, after this, the radical must receive the negative sign, and, 
consequently, x begins to diminish. Again, from the law of continuity, this 
diminution is kept up until x has arrived at the value c. At this point the dimi- 
nution must change into an augmentation, for x cannot fall below c. Thus the 
movement of X is a continuous swinging back and forth between the limits c and b. 

\XT t V (b — x)(x — c) 
We can put X ^^ ^ ^ '- , 

R being a function of jr which remains constantly positive and finite for all values 
of X between c and b. We can then write 

dt R 



dx ^\{b—x){x — c)\ 

A new variable u can now be advantageously introduced in place of x. Let 

X ^ a(\ — e cos li) , 

where a-=-\(b -\- c), and e-=(b — c)\(b -\- c); and u is equivalent to an integral 
number of circumferences when x^c, and augments by half a circumference 
when X, next following, attains the value b. Thus u, like t, augments continu- 
ously. We have 

b — X ■=.ae{\. -|- cos «), x — c ^ ae(i — cos ti) , 
\/^b — x)(x — c)\ = ae sin u, 

dx = ae sin u du. 
Therefore dt = Rdu. 
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As ^ is a one- valued function of ;f or of a ( i — e cos «), it can be expanded 
in the following periodic series 

^ = - [i -f- a, cos u -|- 2«2 cos 2u + 303 cos 3« -f- . . . ], 

n, «,, «2, etc. being constants, the first having the value 



n nj 



o 

Then c being an arbitrary constant, 

« (/ -j- ^) ^ « + a, sin « + Oj sin 2u -f a, sin 3« + . . . . 

This series serves for determining t when x or u is given ; but, more fre- 
quently it is X or « which is required in terms of t. It is necessary, then, to in- 
vert the series. The coefficients of the inverted series are most readily found by 
means of definite integrals. Let us suppose that it is required to find the periodic 
series, in terms of t, for a function of x and Xj which we will denote by U. This 
function we assume to be always finite and continuous. The base of hyperbolic 
logarithms being e, let us put 

fl = n{t-\-c), 2 = £^''-\ s — e'"'-\ 

and, for brevity, 

2S=ai{s — s-^) + a2{s'—s-^) + a3(s^ — s-^) + . . . . 

The equation connecting 2 and s is 



«^ -f- 00 
We can suppose that d^ 2" Ai^'. 



Then A=— CUz-'dr = — Cus-^-^ nRdu. 

2Ttj ^ 2nJ 

o o 

f/ being = F{x, x^, we have 

rT r- i / \ aesinu) 

U ^ F <a{i — ^cos u), — ^ — > 

Supposing that U is reduced to x, it is plain that the coefficient of s^, in the 
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development of x- in powers of z, is the same as the coefficient of s' in the devel- 
opment of 

«[l-i.(5 + .-')][l+.^]s-« 

in powers of s. 

By adopting the Be.sselian functions y^f*', we have 

,-j*..(.-.-.) L I yw),.^ s^, ,-i^,(s'-.-') L I /y),.^ s^, etc.; 
and the expression, given above, can be written 



a < 






+ 



where ^ =: — ^ioj. 

However, unless the coefficients a,, etj, a^, . . . decrease rapidly, this will not 
be a practical method of developing jr in a periodic series. Generally it will be 
shorter to employ mechanical quadratures in obtaining the value of the definite 
integral. Let us suppose that 

•^ = ii^o + A cos C + A cos 2C + A cos 3C + . . . . 
Then 



^i = - \ X cos i!^ d^ 



= - |«(i — £ cos «) [i + «! cos u + 202 cos 2« + 3«3 cos 3« -f . . . ] cos t^ du, 
o 
where, to obtain the value of C corresponding to a given value of u, we employ 
the equation 

ij' z= « -|- Oj sin « + «2 sin 2« + . . . . 

It will be seen that this method is applicable to a much wider range of ques- 
tions than the motion of planets in elliptic orbits. And the superiority of the 
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method of definite integrals over Lagrange's Theorem for the inversion of the 
series is quite manifest. 

II. 

In order to illustrate the preceding general theory, let us treat the problem of 
n material points moving about a centre under the action of central forces admit- 
ting a potential which is a function of the sum of the squares of the radii vectores. 
Each point will then move in a fixed plane and its radius vector will describe 
equal areas in equal times. Thus all will be virtually known in reference to 
these motions, provided we are able to express the radii vectores as functions of 
the time. 

Let the radii be denoted as r,, rj, . . . r„, and the orbit longitudes, measured 
each from any point in its plane, as ^, ^2. • • • ^»- Fo"" brevity, put 

{? = r,' + ri^ . . .+ r„\ 

Then, if the potential is represented by f (/>), we shall have the two equations, 
representing generally all the equations of the problem, 

d^r, ^dK^ _ .r^ 

dk_lu 
dt~ rr 

hi being the constant of areolar velocity. Consequently if we put 

' i 

the general form of the differential equations determining the radii vectores will be 

d\i _ 9ii 

dt' ~ dr, ■ 

They have the integral, corresponding to that of living forces, 

I§=2iii+C), 

C being an arbitrary constant. Also we may derive 

1 r, — r-T = ^ ri X-- . 
* dt'' 3ri 

By adding the last two equations, 



d 
dt 



[^'f]=2f(/')+/'f'(/') + 2C 
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an equation involving only the dependent variable p. Multiplying it by the fac- 
tor 2/> -T^ , and integrating, we get, A being an arbitrary constant. 

Whence ' + ' = i VWmpT^^^Cj^^^^] ' 

Inverting this we shall have /> as a function of t. 

By dividing the penultimate equation by p^ and dtfiferentiating, we get 

dt^~ p '^^ p" ■ 



The general equation determining ^ is 






As p is now a known function of t, r^ is the only unknown in it, and, conse- 
quently, the equation by itself suffices for determining it. To put the equation' 
in a form suitable for integration, let us eliminate {'(p) between the last two 
equations. We get 

d(pdr,-r,dp) _k,^ A' 
dt^ ~ rf^ p" " 

To simplify this, we will adopt an auxiliary variable ^, such that 
'^y-]?d'-p^{2p^[{(p)+(r\-A^}- 



Then 3ii_P^V* ^y^ 

Whence, by integration, we derive 

lMV_2a._tt'-^^^ 
d<p j -^^' p" A'rr 

Bi being the arbitrary constant. By putting 

p •' ^ ' 
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we get dip = 



2^22.,U, — U^ — -^ 



For convenience, adopting a new constant e^, in place of h^, such that hi\JP 
= a<* (i — e^ the quantity under the radical sign becomes 

[a< (i + e^ — «J [«i — aj (i — e^'\ . 

Thus, putting «j = a^ (i — et cos £<), Sj being a new variable, we get d^ = ^dSi, 
and thus s^ = 2^ -|- «<, o^ being a constant. Thus we have, in fine, 

^ = i/{ai [i — ei cos (2<f' + 20^)]}. 

As we have 2* -4 = i , 

P 

the constants a^, i-j, and a, satisfy the relations 

2" a, = I, 2' aj^,- cos 2ai = 0, 2" a<^i sin 2aj = O. 

We thus have 2« independent arbitrary constants introduced by integration ; the 
number there should be. 

In order to find an expression for the longitudes, we take the general equa- 
tion 

hidt 



dh-. 



aip' [i — ei cos 2 {4' + «()] 

i/(i-.,^)4A 



I — Ci cos 2{lp -\- «() ' 

The integral of which gives 

tan (>i, + A) = ^ [f^i] • *^" (^'' + ''^^ ■ 

/9j being the arbitrary constant. 

To simplify the equations which give t + c and tl\ we suppose that a(i + e) 
is the maximum value of p, and a(i — e) its minimum value. Then we can 
adopt a variable e such that 

/> = a ( I — e cos e) . 

Thus df) = a^ sin £</£, and we may put 

2/>' [f (/>) + C] — y42 = /?aV sin^ e, 
where R remains constantly positive throughout the motion of ;0. Then 
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R, being a function of />, is also oneof a(i — f cos s), and thus is capable of being 
expanded in a converging series of terms, each consisting of a constant multi- 
plied by the cosine of a multiple of s. Also (il R and A j pR can be expanded in 
similar series. Then the period T, in which ,« goes through the round of its val- 
ues, is given by the definite integral 



,„ ra ( I — ^ cos s) , 
^=J-^~~R ^^^ 



and the augmentation of the variable (p, in the same time, will be equivalent to 
the definite integral 

2Jr 

Ade 



f: 



a ( I — e cos e) R ' 

u 

If the value of the latter is 2^-, <p will augment by a circumference while ,o goes 
through its period. This is the case when f(/*) = nj p; but, in general, this con- 
dition is not fulfilled. 

Provided that A^ is a positive quantity, it is plain that, after ^ has gone 
through its period, the longitudes and latitudes, whether as seen from the centre 
or from any of the points, all return to the same values. The same thing is true 
of the ratios of the radii vectores. Thus the movement of the system may be 
conceived as taking place under the operation of two distinct causes. The first 
producing a revolution of all the points about the centre in closed curves and in 
the same time, while the second, having a different period, changes the scale of 
representation of the system in space. 

In the preceding treatment we have supposed that A'' is a positive quantity. 
When this is not the case, some modifications must be made. Let us suppose 
first that A^o. Then we have 

dp 



t+c=f. 



l/{2[S{p)+C]}' 

do 



Then J — ^r ( = a^ — ki 



and we may assume </> = J-^^^-^^^£^^_^ 

p-yjl a, J 
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Also dk = ^ = hA^d<p 

hiAidtp 

In the second place, let A^ be negative. Here it is only necessary in some 
places to accomplish the integrations by the aid of hyperbolic cosines instead of 
circular. 

The differential equations of this problem, in the case where the radii are 
supposed to describe no areas, were first integrated by Binet.* But the addition, 
to the forces, of the terms arising from centrifugal action, much enhances the 
interest of the problem. 

*See \AOU\i\\t, Journal Je Mathematiques, First Ser. Tome IL p. 457. 



ON THE FOCAL CHORD OF A PARABOLA. 
By Prof. R. H. Graves, Chapel Hill, N. C. 

Let J/" ^ i,ax be the equation to a parabola, s its focus, and PSP' a focal 
chord. Let the tangent and normal at P' meet the diameter through Pat J/and N. 

It may be easily proved that PM = PN == PP' and that a similar property 
holds for the tangent and normal at P. 

Therefore, if two equal rhombs be constructed on PP' having two other 
sides of each parallel to the axis, their diagonals are tangents and normals at P 
and P' ; and the tangent at one point is parallel to the normal at the other. 

Each normal chord divides the other in the ratio i : 3. 

The chord joining the other ends of the normal chords is parallel to PP' and 
three times as long. 

A line perpendicular to PP' at 5, and terminated by this parallel chord and 
the pole of PP', is divided by 5 in the ratio I : 4. 

Hence the locus of the fool of the perpendicular dropped from 5 on the par- 
allel chord is a right line, whose equation is .r ^ ga. 

Hence the envelope of the parallel chord is a confocal parabola, having for 
its equation /■ = 32^ {ga — x) . 

It cuts the original parabola orthogonally where it is cut by its evolute. 



